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Time-odd mean fields (nuclear magnetism) and their impact on physical observables in rotating 
nuclei are studied in the framework of covariant density functional theory (CDFT). It is shown 
that they have profound effect on the dynamic and kinematic moments of inertia. Particle number, 
configuration and rotational frequency dependences of their impact on the moments of inertia have 
been analysed in a systematic way. Nuclear magnetism can also considerably modify the band 
crossing features such as crossing frequencies and the properties of the kinematic and dynamic 
moments of inertia in the band crossing region. The impact of time-odd mean fields on the moments 
of inertia in the regions away from band crossing only weakly depends on the relativistic mean field 
parametrization, reflecting good localization of the properties of time-odd mean fields in CDFT. 
The moments of inertia of normal-deformed nuclei considerably deviate from the rigid body value. 
On the contrary, superdeformed and hyperdeformed nuclei have the moments of inertia which are 
close to rigid body value. The structure of the currents in rotating frame, their microscopic origin 
and the relations to the moments of inertia have been systematically analysed. The phenomenon 
of signature separation in odd-odd nuclei, induced by time-odd mean fields, has been analysed in 
detail. 
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I. INTRODUCTION 

The development of self-consistent many-body theories 
aiming at the description of low-energy nuclear phenom- 
ena provides the necessary theoretical tools for an ex- 
ploration of the nuclear chart into known and unknown 
regions. Theoretical methods (both relativistic [l[ and 
non-relativistic 0j) formulated within the framework of 
density functional theory (DFT) are the most promis- 
ing tools for the global investigation of the properties 
of atomic nuclei. The power of the DFT models is es- 
sentially unchallenged in medium and heavy mass nuclei 
where 'ab-initio' type few-body calculations are computa- 
tionally impossible and the applicability of the spherical 
shell model is restricted to a few regions in the vicinity 
of doubly shell closures. 

The mean field is a basic concept of every DFT. One 
can specify time-even and time-odd mean fields [H, H| 
dependent on the response of these fields to the action 
of time-reversal operator. While the properties of time- 
even mean fields in nuclear density functionals are rea- 
sonably well understood and defined 0, , there are still 
many unknowns in our knowledge of time-odd mean fields 
which appear only in the nuclear systems with broken 
time-reversal symmetry. This is especially true in the 
covariant density functional theory (CDFT) [l| where 
only few articles were dedicated to the study of time- 
odd mean fields (see Ref. [j| for review). Note that the 
effects, produced by the magnetic potential in the Dirac 
equation and called as nuclear magnetism (NM) [6| in 
the framework of the CDFT, are due to time-odd mean 
fields. 

Rotating nuclei represent a system which is strongly 
affected by time-odd mean fields. The representative 
studies of few examples [1, H], 0-[t| clearly show that the 



kinematic and dynamic moments of inertia of the nuclei 
rotating in collective manner are considerably affected by 
time-odd mean fields. It was shown in the CDFT frame- 
work Q that microscopic mechanism of this modification 
is traced back to the modifications of the expectation val- 
ues of the single-particle angular momentum (j x )i in the 
presence of NM. The contribution to (j x )i due to NM is 
defined as 

a&), = Q x )f M - Qx)Y NM (i) 

where the subscripts NM and WNM indicate the val- 
ues obtained in the calculations with and without NM, 
respectively. The A{j x )i is positive at the bottom and 
negative at the top of the TV-shell [H . The absolute value 
of A(j x )i correlates with the absolute value of (j x )i- Note 
that the contributions to (j x )i due to NM are small in the 
middle of the shell. The A(j x )i contributions can be de- 
composed into the contributions due to spin (A(s x )i) and 
orbital (A(l x )i) angular momenta, which have compli- 
cated dependences both on the frequency and the struc- 
ture of the single-particle orbital under study Q ■ Similar 
features are expected also in non-r elativistic DFT [|. 

The changes in the alignment properties of the single- 
particle orbitals induced by NM (Eq. (Q])) reflect them- 
selves also in physical observables such as effective align- 
ments and the energy splittings between signature part- 
ner orbitals (signature splitting), measured experimen- 
tally Moments of inertia and effective alignments 
in normal- and superdeformed nuclei in different parts 
of nuclear chart [l], l8l-fl4j are well described by the 
parametrizations which include non-linear self-couplings 
only for the tr-meson (see Table 1 in Ref. @]). This fact 
strongly suggests that NM is well accounted in this type 
of the relativistic mean field (RMF) parametrizations. In 
addition, NM can have an impact on the energy gap be- 
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tween the yrast and excited configurations in local min- 
ima (as illustrated on the example of the hyperdeformed 
minima in Ref. 15]), on the terminating states [lH and 
on the additivity of angular momentum alignments [l6j . 

A systematic investigation of time-odd mean fields in 
one- (two-) particle states in odd (odd-odd) non-rotating 
nuclei has been performed in our previous article [j| . The 
current manuscript is a continuation of our efforts aimed 
on comprehensive understanding of time-odd mean fields 
in CDFT. Its goal is a systematic study of time-odd mean 
fields and their manifestation in rotating nuclei. 

Tabic 1 in Ref. @ shows large variety of the 
parametrizations of the RMF Lagrangian. The investi- 
gation of all these parametrizations is definitely beyond 
the scope of this study. Thus, the present investigation 
has been focused on the study of time-odd mean fields 
in the CDFT with the parametrizations of the RMF La- 
grangian including only non-linear self-couplings of the 
cr-meson (the group A of the parametrizations in Table 
1 of Ref. [5j). So far, only this group of parametriza- 
tions has been used in the study of rotating nuclei [l|, ®- 
UH E3, EH ■ The results of the study of time-odd mean 
fields in the groups B, C, and D of the parametrizations 
of meson-coupling models (Table 1 in Ref. [Hj]) as well 
as within the point-coupling models will be presented in 
a forthcoming manuscript. 

The manuscript is organized as follows. The cranked 
RMF theory and its details related to time-odd mean 
fields in rotating nuclei are discussed in Sec. [TXJ Section 
IIIII is devoted to the analysis of the impact of time-odd 
mean fields on band crossing features. Particle number 
and deformation dependences of the impact of NM on the 
moments of inertia are considered in Sec. IIVI Currents 
(and their single-particle origin) in the intrinsic frame of 
rotating nuclei are discussed in Sec. [V] Frequency and 
configuration dependences of the impact of NM on the 
moments of inertia are analyzed in Sec. lVIl Parametriza- 
tion dependence of the NM contributions to the moments 
of inertia are discussed in Sec. lVIII Sec. lVIIIl is devoted to 
the study of time-odd mean fields in terminating states. 
The phenomenon of signature separation in odd-odd nu- 
clei is investigated in Sec. lIXI Finally, Sec. |X] reports the 
main conclusions of our work. 



II. THEORETICAL FORMALISM 

The results presented in the current manuscript have 
been obtained mainly in the framework of Cranked Rel- 
ativistic Mean Field (CRMF) theory [|, H ||. Only 
in a few cases the results obtained within the Cranked 
Relativistic Hartree-Bogoliubov (CRHB) theory [l3j are 
shown. The CRMF theory has been successfully em- 
ployed for the description of rotating nuclei (see Refs. 
[l|, [l8| and references therein) in which time-odd mean 
fields play an important role. In this theory the pairing 
correlations are neglected which allows to better isolate 
the effects induced by time-odd mean fields. The most 



important features of the CRMF formalism related to 
time-odd mean fields are outlined below (for more details 
see Refs. [1,0), while the details of the CRHB theory are 
presented in Ref. (l3| . 

In the Hartree approximation, the stationary Dirac 
equation for the nucleons in the rotating frame (in one- 
dimensional cranking approximation) is given by 



(h D - VL X J X )'4> 1 = Eiipi 



(2) 



where hp is the Dirac Hamiltonian for the nucleon with 
mass to 

h D =a(-iV -V(r)) + Vo(r) + p(m + S(r)) (3) 
and the term 
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is just the Coriolis term. Note that the rotational fre- 
quency Q x along the s-axis is defined from the condition 
that the expectation value of the total angular momen- 
tum at spin I has a definite value 



J(Sl x ) = <$ n I Jx I $n> = y/ 1(1 + 1)- 



(5) 



The Dirac Hamiltonian contains the average fields de- 
termined by the mesons, i.e. the attractive scalar field 
S(r) 



S(r) = ga<j(r), 



(6) 



and the repulsive time-like component of the vector field 
Voir) 



V (r) = g u uj a (r) + g p T 3 p {r) + e 
A magnetic potential V(r) 

V(r) = g u u>(r) + g p T 3 p(r) 



1-T 3 



Mr)- (7) 



A(r) 



(8) 



originates from the space-like components of the vec- 
tor mesons. Note that in these equations, the four- 
vector components of the vector fields w M , p^, and 
are separated into the time- like (ujq, po and Aq) and 
space-like [u> = (oj x , uj v , uj z ) , p = (p x , p v , p z ), and A = 
(A x ,A y ,A z )] components. In the Dirac equation the 
magnetic potential has the structure of a magnetic field. 

The corresponding meson fields and the electromag- 
netic potential are determined by the Klein-Gordon equa- 
tions 

{-A + to^} a (r) = -g a [p n s (r) + pl{r)] 

-g2<J 2 (r)~g 3 a 3 (r), (9) 
{~A + ml}oj (r) = 9 M(r) + pP(r)}, (10) 
{-A + ml} u (r) = g^j n (r)+j p (r)} (11) 
|-A + TO2} Po (r) = g p [p n v {r) - p%{r)l (12) 
{-A + to^} p(r) = 9 P [j n {r)-f{r)l (13) 
-AAo(r) = epP(r), —AA(r) = ef(r), (14) 
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with source terms involving the various nucleonic densi- 
ties and currents 

N,Z 

P^ p (r) = £>i(r))Wi(r), (15) 
t=i 

N,Z 

P n v ' P (r) = £(^(r))t^(r), (16) 

i=l 
AT,.Z 

r p (r) = j2(Mr)?&Mr) (17) 

i=l 

where the labels n and p are used for neutrons and pro- 
tons, respectively. In the equations above, the sums run 
over the occupied positive-energy shell model states only 
(no-sea approximation) [20l . |2l| . Note that the spatial 
components of the vector potential A(r) are neglected in 
the calculations since the coupling constant of the electro- 
magnetic interaction is small compared with the coupling 
constants of the meson fields. 

Two terms in the Dirac equation, namely, the Corio- 
lis operator J x and the magnetic potential V(r) (as well 
as the currents j n ' p {r) in the Klein-Gordon equations) 
break time-reversal symmetry Their presence leads 
to the appearance of time-odd mean fields. However, one 
should distinguish time-odd mean fields originating from 
Coriolis operator and magnetic potential. The Coriolis 
operator is always present in the description of rotating 
nuclei in the framework of the cranking model. However, 
the CRMF calculations, with only these time-odd fields 
accounted for, underestimate the experimental moments 
of inertia [8, 9]. A similar situation also holds in non- 
relativistic theories [!, [22| ■ The inclusion of the currents 
j n ' p (r) into the Klein-Gordon equations, which leads to 
the space-like components of the vector w and p mesons 
and thus to magnetic potential V(r), considerably im- 
proves the description of experimental moments of iner- 
tia. The effect coming from the space-like components 
of the vector mesons is commonly referred to as nuclear 
magnetism @ since the magnetic potential has the struc- 
ture of a magnetic field in the Dirac equation. 

Note that time-odd mean fields related to NM are de- 
fined through the Lorentz invariance [l[ and thus they 
do not require additional coupling constants: the cou- 
pling constants of time-even mean fields are used also for 
time-odd mean fields. 

The goal of the current manuscript is to understand the 
impact of nuclear magnetism (NM) on the properties of 
rotating nuclei. We will use the terms nuclear magnetism 
and time-odd mean fields interchangeably throughout the 
manuscript. However, one should keep in mind that the 
latter term is related only to the time-odd mean fields 
produced by the magnetic potential. 

Single-particle orbitals are labeled at rotational fre- 
quency fl x = 0.0 MeV by [Nn z A]fl sl9n . [Nn z A]fl arc the 
asymptotic quantum numbers (Nilsson quantum num- 
bers) of the dominant component of the wave function. 
The superscripts sign to the orbital labels are used some- 



times to indicate the sign of the signature r for that or- 
bital (r = ±i). Note that the labelling by means of Nils- 
son labels is performed only when the calculated shape 
of nuclear configuration is prolate or near-prolate. 

Many-particle configurations (further nuclear configu- 
rations or configurations) are specified by the occupation 
of available single-particle orbitals. In the calculations 
without pairing, the occupation numbers n are integer 
(n = or 1). In addition, in the CRMF code it is possi- 
ble to specify the occupation of cither r — +i or r = —i 
signature of the single-particle state. In odd-odd nuclei, 
all single-particle states of specific (proton and neutron) 
subsystem with exception of one are pairwise occupied. 
We will call this occupied single-particle state of fixed 
signature for which its time-reversal (signature) counter- 
part state is empty as blocked state in order to simplify 
the discussion. The specification of nuclear configuration 
by means of listing all occupied single-particle states is 
unpractical. Thus, in odd-odd nuclei the Nilsson labels 
of the blocked proton and neutron states and their sig- 
natures are used for configuration labelling. 

The CRMF equations are solved in the basis of 
an anisotropic three-dimensional harmonic oscillator in 
Cartesian coordinates characterized by the deformation 
parameters /3o and 7 as well as the oscillator frequency 
fajQ = 41 A -1 / 3 MeV. Our selection of the deformation 
parameters of the basis is the same as in earlier system- 
atic studies of rotating nuclei in different regions of nu- 
clear chart (Refs. (9Hlll[T3,[ll). 7 = 0° is used in all cal- 
culations. The deformation parameter /3q of the basis is 
selected in such a way that it provides the convergence to 
the local minimum under study. Thus, /3q = 0.25 is used 
in the case of normal-deformed (ND) states, (3q = 0.2 
and Po = 0.5 in the case of superdeformed (SD) states 
in the A ~ 60 and A ~ 150, 190 mass regions, respec- 
tively, and (3q — 1.0 in the case of the hyperdeformed 
(HD) states. The truncation of basis is performed in 
such a way that all states belonging to the shells up 
to fermionic Np=12 and bosonic Nb=20 are taken into 
account in the calculations of ND and SD states. The 
fermionic basis is increased up to Np—14 in the calcu- 
lations of HD states. Numerical analysis indicates that 
these truncation schemes provide sufficient numerical ac- 
curacy for the physical quantities of interest (see Refs. 

The majority of the calculations are performed with 
the NL1 parametrization [2l[ of the RMF Lagrangian 
since this parametrization has been used extensively in 
earlier systematic studies of rotating nuclei across the 
nuclear chart [jJO, B Q El • 

To investigate the impact of NM (time-odd mean 
fields) on physical observables, the CRMF calculations 
are performed in two calculational schemes for fixed con- 
figurations: 

• Fully self-consistent calculations with NM included 
(hereafter denoted NM calculations), which take 
into account space-like components of the vector 
mesons (Eqs. (TTTj) . (|T3"|) and flSJ)), currents (Eqs. 
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(fTTj) . (fT3)l . and ([T7))). and magnetic potential V(r) 
(Eq. ©). 

• Fully self-consistent calculations without NM 
(hereafter denoted as WNM calculations), which 
omit space-like components of the vector mesons 
(Eqs. ([HI), (US]) and ©), currents in the Klein- 
Gordon equations (Eqs. (fTTj) and (fT5|) . and mag- 
netic potential V(r) (Eq. (JSj) ) . The results of the 
NM and WNM calculations are always compared 
for the same nuclear configuration. 

These are the ways in which the effects of time-odd 
mean fields can be studied, and as such they are fre- 
quently used in DFT studies of rotating systems, both in 
relativistic and non-relativistic frameworks [3, 0, & EH, 
|22T |. One should, however, keep in mind that if time- 
odd fields are neglected, the local Lorentz invariance 
(Galilean invariance in non-relativistic framework 0,[23|) 
is violated. The inclusion of time-odd mean fields restores 
the Lorentz invariance. 

It is interesting to compare the basic features such as 
Lorentz invariance and the definition of the coupling con- 
stants of the time-odd channel of the CDF theory dis- 
cussed above with the ones of non-relativistic Skyrme en- 
ergy density functional (EDF) theory. It was recognized 
in earlier Skyrme DFT studies, that the connection be- 
tween the coupling constants of time-odd and time-even 
channels depends on what entity, namely, Skyrme force 
or energy density functional is considered to be more fun- 
damental 

H SI Hi- H the Skyrme force is considered 
more fundamental then the time-odd constants are de- 
termined as a function of time-even constants [H, l24j . 
However, since the time-even coupling constants are usu- 
ally adjusted solely to the time-even observables, the re- 
sulting values of the time-odd coupling constants simply 
"fictitious" or "illusory", as noted already in Ref. [25| 
(see also Ref. [13] )■ On the contrary, in the framework 
of the Skyrme energy density functional theory, time-odd 
properties of the functional are independent of time-even 
properties which is a consequence of broken link between 
the Skyrme force and the density functional. 

The question of whether Galliean invariance must be 
imposed in Skyrme EDF is not yet resolved (23^ . despite 
the fact that it is imposed in many studies. Note that in 
many phenomenological approaches, such as the nonin- 
teracting or interacting shell models, Galilean symmetry 
is not considered, because the translational motion is not 
within the scope of such models [23[ . 

It is also important to mention that the cranking mod- 
els based on phenomelogical Woods-Saxon or Nilsson po- 
tentials do not incorporate time-odd mean fields. How- 
ever, they succesfully describe rotating nuclei [26l W\. 

Note that the Coriolis term is present in NM and WNM 
calculations. This means that the currents (Eq. (|T7|) ) are 
always present in rotating nuclei. However, it is impor- 
tant to distinguish the currents induced by the Coriolis 
term and the ones which appear due to magnetic po- 
tential. The currents, which appear in the WNM cal- 



culations, are generated by the Coriolis term. Thus, we 
will call them as Coriolis induced currents. On the con- 
trary, the currents in the NM calculations are generated 
by both the Coriolis term and magnetic potential. The 
difference of the currents in the NM and WNM calcu- 
lations is attributable to magnetic potential. Thus, the 
currents \j n ' v (r)] NM - [j n ' p (r)] WNM will be called mag- 
netic potential induced currents. 

In the following, the contribution AO NM ~ contr (in per- 
centage) of NM to the physical observable O is defined 
as 



AO 



NM — contr 



o 



NM 



o 



WNM 



o 



NM 



x 100%. 



(18) 



The physical observables, most frequently used in the 
analysis of rotating nuclei, are kinematic (J^) and dy- 
namic (J^ 2 ') moments of inertia which are defined as 



J w (n x ) = ±, 



where J is the expectation value of the total angular 
momentum along the x-axis. In the CRMF theory, this 
quantity is defined as a sum of the expectation values of 
the single-particle angular momentum operators j x of the 
occupied states 



(20) 



Thus, the modifications of the moments of inertia due 
to NM can be traced back to the changes of the single- 
particle expectation values (j x )i = (i\]x\i) and the corre- 
sponding contributions of spin ((s x )i) and orbital ({l x )i) 
angular momenta 



III. NUCLEAR MAGNETISM AND BAND 
CROSSING FEATURES 

Since NM substantially modifies the single-particle 
properties (energies, alignments) [3, 0|, it is reasonable 
to expect that the band crossing features are affected by 
NM. In order to study this question, the CRMF (without 
pairing) and the CRHB+LN calculations have been per- 
formed for lowest superdeformed (SD) band in 194 Pb. In 
the CRHB+LN calculations, the D1S Gogny force [H] 
is used in pairing channel and an approximate particle 
number projection is performed by means of the Lipkin- 
Nogami method (T3 |. 

The unpaired proton band crossing seen in the CRMF 
calculations originates from the interaction between the 
tt[642]5/2+ and tt[651]1/2+ orbitals (Fig.UJi). Since NM 
increases somewhat the single-particle alignment (j x )i 
(Fig. [lb) and the slope of the routhian for the 7r[651]l/2+ 
orbital (Fig. [TJi), the band crossing takes place at lower 
frequency. The shift of crossing frequency due to NM is 
considerable (120 keV) from 0.465 MeV (WNM) down to 
0.345 MeV (NM), Fig. [TJi. The calculations also suggest 
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NM WNM 




Rotational frequency £2 X [MeV] 



FIG. 1: (a) Proton single-particle energies (Routhians) in the 
self-consistent rotating potential as a function of rotational 
frequency Q x obtained in the CRMF calculations with and 
without NM. They are given along the deformation path of the 
lowest SD configuration in 194 Pb. Only interacting [651]l/2 + 
and [642]5/2 + orbitals are shown, see Fig. 1 in Ref. for 
full spectra, (b) The expectation values (jx)% of the single- 
particle angular momentum operator j x of the orbitals shown 
on panel (a). Solid and dashed arrows are used to indicate the 
frequencies (as well as the energy gap between the interacting 
orbitals in panel (a)) at which the band crossings take place 
in the calculations with and without NM, respectively. 
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FIG. 2: Kinematic (J^) and dynamic (J 1 - 2 ') moments of 
inertia for the lowest SD configuration in 194 Pb obtained in 
the calculations with and without NM. Proton and neutron 
contributions to these quantities are indicated by ir and v, 
while total moments by Hot'. Panels (a) and (b) show the re- 
sults obtained in the calculations without pairing, while pan- 
els (c) and (d) show the results of the calculations within the 
CRHB+LN framework. Solid and dashed arrows are used 
to indicate the frequencies at which the band crossings take 
place in the calculations with and without NM, respectively. 



that the strength of the interaction between two inter- 
acting orbitals at the band crossing is modified in the 
presence of NM as seen in the change of the energy dis- 
tance (gap) between these two orbitals at the crossing 
frequency (Fig. UK). 

An additional mechanism affecting the band crossing 



frequencies will be active in odd- and odd-odd mass nu- 
clei as well as in excited configurations of even-even nu- 
clei. In such configurations, there is at least one single- 
particle state the opposite signature of which is not oc- 
cupied. This results in the currents at Vl x = 0.0 MeV 
,5:j. The energy splitting between different signatures of 
the single-particle states at no rotation is a typical conse- 
quence of these currents (see Sec. IVA in Ref. [5| for more 
details). As a result, the energy gap between interacting 
orbitals at D, x = 0.0 MeV can become larger or smaller 
dependent on the impact of the currents on the single- 
particle energies of interacting states. Consequently, this 
change in the energy gap will translate into higher or 
lower band crossing frequencies. Note that for simplicity 
we assume that the Q x = 0.0 currents will not modify the 
alignment properties of interacting orbitals; this trans- 
lates into the independence of the single-particle routhian 
slope in the energy versus Q x plot (see, for example, Fig. 
QJ,) on the Q x = 0.0 currents. 

Fig. UK can be used to illustrate this mechanism. Let 
assume that the fl x = 0.0 currents will increase the en- 
ergy gap between the 7r[642]5/2 + and 7r[651]l/2 + orbitals 
at Q x = 0.0 MeV: this will lead to higher band crossing 
frequencies. However, the band crossing frequencies will 
decrease in the case when the energy gap between these 
orbitals at Q x = 0.0 MeV becomes smaller in the pres- 
ence of the Q x — 0.0 currents. The assumption that the 
£l x = 0.0 currents do not have an impact on the align- 
ment properties of interacting orbitals is definitely too 
simplistic but it allows to illustrate the fact that NM can 
both decrease and increase the band crossing frequencies. 

This mechanism is not active in the configuration of 
even-even 194 Pb nucleus discussed above since both sig- 
natures of all states below the Fermi level are pairwise 
occupied. As a result, no current is present at Q x = 0.0 
MeV. 

The impact of NM on band crossing features is also 
seen in the CRHB+LN calculations where the alignment 
of the pairs of J15/2 neutrons and protons causes the 
shoulder and peak in total dy namic moment of inertia 
J( 2 ) (Fig. [5b) (see also Ref. [Taj), Note that each of these 
two alignments creates a peak in the dynamic moment 
of inertia of corresponding subsystem. NM shifts the 
paired neutron band crossing to lower frequencies by 70 
keV from 0.485 MeV (WNM) to 0.415 MeV (NM). Paired 
proton band crossing lies in the calculations with NM 
at fl x = 0.535 MeV, while only the beginning of this 
crossing is seen in the calculations without NM (Fig. [2b). 

The origin of this effect is twofold. Similar to the un- 
paired calculations, the part of it can be traced to the 
fact that NM increases the expectation values (j x )i of 
the orbitals located at the bottom of the shell (the dis- 
cussed orbitals are of this kind) [J|. The corresponding 
larger slope of the quasiparticle routhians causes the shift 
of the crossing to lower frequencies. However, an addi- 
tional contribution comes from the modification of the 
pairing by NM. There is a difference in the pairing en- 
ergies calculated with and without NM which increases 
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FIG. 3: Transition quadrupole moments Q t (panel (a)), mass 
hexadecapole moments Q40 (panel (b)), proton and neutron 
pairing energies E pa iring = — 1/2 Tr(An) (panel (c)) and pro- 
ton and neutron particle number fluctuations {(AN) 2 ) (panel 
(d)) of the lowest SD solution in 194 Pb obtained in the calcu- 
lations with and without NM. Thick and thin lines in the up- 
per panels are used for the results obtained in the CRHB and 
CRMF theories, respectively. On the bottom panels, thick 
and thin lines are used for neutron and proton quantities, 
respectively. 



with rotational frequency, see Fig. [3]:,d. The pairing in 
the calculations with NM is weaker. This can be ex- 
plained by the increase of (j x )i of the orbitals located 
at the bottom of the shell due to NM (see above). The 
gradual breaking of high-j pairs proceeds faster, which 
is reflected in a faster decrease of pairing with increasing 
Q x . Thus we can specify this effect as an anti-pairing 
effect induced by NM. 

These considerable differences in the crossing frequen- 
cies obtained in the calculations with and without NM 
cannot be attributed to the differences in equilibrium de- 
formations, since calculated transition quadrupole mo- 
ments Qt and mass hexadecapole moments Q40 obtained 
in the calculations with and without NM differ only 
marginally before band crossing, see Figs. [3k,b. 

The influence of time-odd mean fields on band crossing 
features has been studied by means of a schematic non- 
self-consistent model based on the Skyrme forces in Ref. 
0. In this study, time-odd fields emerging from the S 2 
and — SAS terms of the Skyrme Hamiltonian shift the 
alignment of the 113/2 neutron pair to higher frequen- 



cies in 5S Dy. On the contrary, this crossing appears at 
lower frequencies in the CRHB+LN calculations when 
NM is taken into account. This difference is not surpris- 
ing considering the fact that time-odd mean fields are not 
well defined in non-relativistic density functional theo- 
ries [3, H^] . It was also suggested in the cranked Skyrme 
Hartree-Fock framework that time-odd mean fields may 
be responsible for band crossing in yrast superdeformed 
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FIG. 4: (Color online) The contribution (in %) of NM to 
the kinematic moments of inertia of nuclei in different isotope 
chains at normal deformation. The results for the lowest in 
energy solutions are shown at rotational frequency Q x — 0.3 
MeV. The frequency has been fixed at this value in order to 
make the comparison with the results of Ref. straightfor- 
ward. 



band of 60 Zn [30]. However, this crossing is described as 
paired band crossing in the CRHB+LN calculations [3lj |. 

Above discussed CRMF and CRHB+LN examples 
clearly show that the modifications of band crossing fea- 
tures (crossing frequencies and the features of the kine- 
matic and dynamic moments of inertia in band cross- 
ing region) caused by NM are substantial and depend on 
the underlying modifications of single-particle properties 
such as alignments and single-particle (quasi-particle) en- 
ergies. 



IV. PARTICLE NUMBER AND DEFORMATION 
DEPENDENCES OF THE IMPACT OF 
NUCLEAR MAGNETISM ON THE MOMENTS 
OF INERTIA 

In the current section, the particle number and defor- 
mation dependence of the impact of NM on the kinematic 
moments of inertia are discussed in detail. We consider 
the contribution of NM to kinematic moment of inertia, 



namely, the (J 



(i) 

NM ' 



' j wnm)I j nm quantity, and its vari- 
ations as a function of particle number and deformation. 
In addition, we investigate how close fully self-consistent 
value of the kinematic moment of inertia comes to the 
rigid body moment of inertia J r i g - The latter quantity 
is obtained in one-dimensional cranking approximation 
with the rotation defined around the x-axis from the cal- 
culated density distribution p(r) by 



p(v)(y 2 + z 2 )d 3 r 



(21) 



The contributions of NM to kinematic moment of in- 
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ertia (the (J$' M - Jwnm)I J nm quantity) for normal 
deformed bands in a number of isotope chains with pro- 
ton number Z > 50 are shown as a function of neutron 
number in Fig. 2J Only the cases in which the nuclear 
configurations are the same in the calculations with and 
without nuclear magnetism are shown in this figure. NM 
typically increases the calculated kinematic moments of 
inertia by 10-30%. However, this increase is around 40% 
in the N — 108, 110 W isotopes. Considerable fluctua- 
tions of the (Jjjm ~ Jwnm)/Jnm quantity as a func- 
tion of neutron number seen in some isotope chains are 
due to the changes in underlying single-particle struc- 
ture. Large changes in the (J]p M - J$nm)/ J nm quan- 
tity are seen on going from the isotope with neutron 
number N to the isotope with N + 2 when two neu- 
tron single-particle orbitals, by which the configurations 
of compared nuclei differ, have the expectation values 
of the single-particle angular momenta (j x )i strongly af- 
fected by NM. The opposite is also true when two neu- 
tron single-particle orbitals, by which the configurations 
of compared nuclei differ, have the expectation values of 
the single-particle angular momenta (j x )i that are only 
marginally affected by NM. Note that in some cases pro- 
ton configurations of two neighboring nuclei with neutron 
numbers N and N + 2 are also different due to the de- 
formation changes; this also contributes into the fluctua- 
tions of the (J?P M — Jwnm)/Jnm quantity as a function 
of neutron number. 

One can also extract from Fig. [4] the dependence of 
the contributions of NM to kinematic moments of in- 
ertia on proton number Z by considering the results 
of the calculations at constant value of neutron num- 
ber N. Such analysis reveals the fluctuations in the 
(Jnm ~ J\vnm) I ^NM quantities which are similar to the 
ones discussed above. The origin of these fluctuations 
can again be traced back to the changes (as a function of 
proton number) in underlying single-particle structure. 

Fig. [5] compares rigid body moments of inertia J r i g 
(Eq. (|2ip) with fully microscopic kinematic moments of 
inertia J^ M (Eqs. (|T9|) and (|2"U|) ) obtained in the calcula- 
tions with NM using the ( J r i g — J^\j)/ Jrig quantity. One 
can see that considerable deviations (in majority of the 
cases being in the window of ±30% but reaching ±60% 
in some nuclei) between these two moments of inertia are 
observed at normal deformation. 

The analysis within the framework of the periodic orbit 
theory [32] concluded that the deviations of the moments 
of inertia from the rigid-body value at high spin are de- 
termined by the shell structure of a system of indepen- 
dent fermions confined by a leptodermous potential. For 
the case of prolate deformation and the rotation perpen- 
dicular to the symmetry axis (the majority of the cases 
studied in the current manuscript fall under this cate- 
gory), the meridian orbits determine the shell moments 
of inertia because only they enclose rotational flux [32| • 

Large similarities are seen between the results of our 
calculations and the ones based on the cranked Woods- 
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FIG. 5: (Color online) The difference (in %) between the 
rigid body moments of inertia J r ; 9 and kinematic moments 
of inertia calculated with NM for the nuclear configurations 
shown in Fig. f3] 



Saxon potential in Ref. |32j ]. 
torn panel in Fig. 10 of Ref. 



Fo r example, right bot- 
shows the difference 
calcu- 



32] 



Jpper — Jrig between the moments of inertia Jp per 
lated in the cranked Woods-Saxon potential and rigid- 
body moments of inertia J r i g for the case of prolate 
deformation and the rotation around the axis perpen- 
dicular to the symmetry axis. If one corrects for the 
difference in the representation of calculated quantities 

[{Jpper - Jrig) in Ref. [HJ and (J ng - jij M )/J r ig in the 
present manuscript], then one can see that our results 

show similar shell dependence of the ( jjyL — J^ig) quan- 
tities as the one seen in Fig. 10 of Ref. (32[. Some dif- 
ferences between these two calculations are in part due 
to simplistic method of the calculation of the rigid-body 
moments of inertia in Ref. [HJ (see Sec. IIB of Ref. [32[ 
for details). 

The CRMF calculations describe rather well the kine- 
matic moments of inertia of normal-deformed [HI, Hl[ 
and smooth-terminating [l], [33| bands at high spin where 
the pairing is negligible. Experimental data on kine- 
matic moments of inertia of normal-deformed rotational 
bands at low spin [which are strongly affected by pair- 
ing] are also well described in the cranked relativis- 
tic Hartree-Bogoliubov calculations [lj], HiJ. These re- 
sults together with the ones presented in the current 
manuscript strongly support the conclusion that weakly- 
and normal-deformed nuclei show the moments of iner- 
tia which strongly deviate from the rigid-body value (see 
also Refs. [SlUHf). 

Figs. [5] and [7] show the results of calculations for yrast 
SD configurations in the A ~ 150 mass region of superde- 
formation and for yrast hyperdeformed (HD) configura- 
tions in the Z = 40 — 58 part of the nuclear chart, respec- 



tively. It is clearly seen that the (J 
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(Color online) The contribution (J NM — 
) (in %) of NM to the kinematic moments of 



inertia and the difference (J r i g — J nm) I ^ ri a ( m %) between 
the rigid body moments of inertia and kinematic moments of 
inertia calculated with inclusion of NM. The results for yrast 
SD configurations in the A ~ 150 mass region of superdefor- 
mation are shown at rotational frequency Q. x = 0.6 MeV for 
nuclei which werepreviously analyzed in the CRMF calcula- 
tions in Refs. ITol]. 




60 

Neutron number N 

FIG. 7: (Color online) The same as in Fig. [6] but for the 
hyperdeformed configurations. The results of the calculations 
are shown for the yrast HD configurations, studied in the 
systematic survey of the hyperdeformation in the Z = 40 — 58 
region of the nuclear chart [TH . [3], at spins at which they 
become yrast (rotational frequency Q, x ~ 0.8 — 1.0 MeV). 



the current manuscript. 



and (J r ig — J nm) I Jrig quantities at these extreme de- 
formations show much smaller fluctuations than the ones 
at normal deformation. Indeed, the contribution of NM 
into kinematic moment of inertia at SD and HD is in 
narrow 20 — 27% range (Figs. [6] and [7]), while it covers 
much large 9 — 43% range at normal deformation (Fig. 
2]). In addition, the values of kinematic moment of in- 
ertia calculated with NM are typically within 5% of the 
rigid body value for the moment of inertia at SD and HD 
(Figs. [6] and [7]), while much larger fluctuations (typically 
within 40% of the rigid body value) are seen in the case 
of normal deformation (Fig. [SJ. 

Microscopic origin of these features can be traced back 
to the underlying shell structure. The analysis within the 
periodic orbit theory [32[ shows that the single-particle 
orbits that cause shell structure of prolate superdeformed 
nuclei do not carry rotational flux if the axis of rotation 
is perpendicular to the symmetry axis. Therefore, the 
moments of inertia of the SD bands in such nuclei should 
be equal to the rigid body value [HJ. Such conclusion 
is in general supported by our microscopic calculations 
which show that the calculated moments of inertia are 
typically within 5% of rigid-body value. The experimen- 
tal deviations (obtained under spin assignments of Refs. 
[HI HII) from the rigid-body values are about 6% or less 
in the A ~ 150 region of superdeformation (see Fig. 5 in 
Ref. [H). 

We also expect that similar mechanism is also responsi- 
ble for the observed features of the moments of inertia at 
HD. However, the periodic orbit theory analysis of such 
features is not available and it goes beyond the scope of 



V. CURRENTS IN INTRINSIC (ROTATING) 
FRAME OF COLLECTIVELY ROTATING 
NUCLEI 

Current distributions in the intrinsic (rotating) frame 
have been studied earlier in several publications. It is well 
known that there are no currents in the intrinsic frame if 
the rigid non-spherical body rotates uniformly (rigid ro- 
tation) (see Sec. 6A-5 in Ref. [36]). The general aspects of 
the velocity (current) fields have been discussed in detail 
in the framework of single-particle Schrodinger fluid [37| , 
which exhibits a remarkably rich variety of fluid dynami- 
cal features, including compressible flow and line vortices. 
Nuclear intrinsic vorticity and its coupling to global rota- 
tions have been studied within the so-called routhian ap- 
proach both in semiclassical approach (38l.l39j and in fully 
self-consistent cranked Hartree-Fock and Hartree-Fock- 
Bogoliubov approaches based on the Skyrme force [4fj| . 
The current distributions in rotating frame have been 
studied in phenomenological cranking approaches based 
on harmonic oscillator [H, |41H431 | and Nilsson [44[ poten- 
tials and in self-consistent cranking approaches based on 
the Skyrme force [40|, |45j . Note that the intrinsic current 
field (as any vector field according to the Hemholtz's the- 
orem) can be split into irrotational and intrinsic vortical 
fields H|. 

Fig. [3] shows typical current distributions obtained in 
the CRMF calculations for selected normal-, super- and 
hyperdeformed nuclei. Despite the fact that the moments 
of inertia of the SD and HD configurations are very close 
to the rigid-body values (Sec. IIV[) . the presence of strong 
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FIG. 8: (Color online) Total neutron current distributions j n (r) in the intrinsic frame in the y — z plane for several normal- 
deformed (ND) (upper row), superdeformed (SD) (middle row) and hyperdeformed (HD) (bottom row) configurations in 
different nuclei. They are shown at x w 0.48 fm in the case of ND and SD configurations, and at x ~ 0.42 fm in the case of 
HD configurations. The results of calculations are shown at rotational frequencies Q. x — 0.3 MeV and Sl a — 0.5 MeV for the 
ND and SD configurations, respectively, and at the spin (Q x ~ 1.0 MeV) at which the HD configurations become yra.st (see 
Ref. pi| for details) in the case of HD configurations. The currents in panels (d-g) and (i) are plotted at arbitrary units for 
better visualization. The currents in other panels are normalized to the currents in panels (d-g) and (i) by using factor F. This 
factor is chosen in such way that the current distribution for every nucleus is clearly seen. The shape and size of the nucleus 
are indicated by density lines which are plotted in the range 0.01 — 0.06 fm -3 in step of 0.01 fm -3 . 



vortices 1 demonstrates the dramatic deviation of the cur- 
rents from rigid rotation. For example, the HD configu- 
rations in 92 Mo (Fig. [Sk.) and 108 Cd (Fig. [SJd) show two 
strong vortices centered at z « ±2 fm. Note that the vor- 
tices (i.e. the curl) of the current fields are aligned or an- 
tialigned along a principal x-axis of the ellipsoid because 
of the use of one-dimensional cranking approximation. 
On the other hand, the HD configuration in 118 Te (Fig. 
[8}:) shows one very strong vortice centered at z = fm, 
and 2 weaker vortices centered atzw ±4.5 fm. All three 
vortices rotate clockwise. The currents in the rotating 
frame of reference that is fixed to the body are caused by 
quantized motion of the fermions. Thus, the differences 
between the currents in 92 Mo and 108 Cd on one hand and 



The existence of vortices at these points implies non-vanishing 
current circulations which are defined as C(r) = V X j(r) [43l |. 



the ones in 118 Te on the other hand are caused by the dif- 
ferences in the underlying single-particle configurations. 
Contrary to the HD configurations, the current distribu- 
tions in the SD configurations of 142 Sm, 148 Gd and 152 Dy 
are characterized by a single very strong central vortice 
(Figs.[5Ji-f). Current patterns in normal deformed nuclei 
100 Sn and 118 Ba look more disordered than in the SD and 
HD nuclei (Figs.[5g,h). This is because three (four) large 
vortices in 100 Sn ( 118 Ba) are spread out over the volume 
of the nucleus. On the other hand, the current pattern 
is dominated by a single large central vortice in the ND 
configuration of 136 Nd (Fig.0). 

Note that all considered configurations are character- 
ized by the weak current in the surface area. On the 
contrary, the average intrinsic current flows mainly in the 
nuclear surface in the semiclassical description of currents 
in normal and superfluid rotating nuclei [38| ■ This under- 
lines the importance of quantum mechanical treatment of 
the currents. 
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FIG. 9: (Color online) Current distributions j n (r) produced by single neutron in indicated single-particle states of the yrast 
SD configuration in 152 Dy at rotational frequencies Q x — 0.0 MeV (left panels) and Q, w = 0.5 MeV (right panels). The shape 
and size of the nucleus are indicated by density line which is plotted at p — 0.01 fm -3 . The currents in panels (e),(f) are 
plotted at arbitrary units for better visualization. The currents in other panels are normalized to the currents in panels (e) 
and (f) by using factor F. The currents and densities are shown in the intrinsic frame in the y — z plane at x = 0.48 fm. The 
single-neutron density distributions due to the occupation of the indicated Nilsson state are shown by colormap. Note that 
slightly different colormap is used in panel (a) for better visualization of densities. 
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The total current is the sum of Coriolis induced and 
magnetic potential induced currents. Total current is 
dominated by the Coriolis induced currents; magnetic po- 
tential induced currents represent approximately 5-20% 
[30%] of total current in the HD and SD [ND] nuclei 
shown in Fig. HI The only exception is 92 Mo, in central 
region of which magnetic potential induced currents are 
larger than Coriolis induced currents by a factor close to 
2. The spatial distribution of Coriolis induced and mag- 
netic potential induced currents is similar in the majority 
of nuclei shown in Fig. [5] However, the spatial distribu- 
tion of these two types of currents differ substantially in 
92 Mo, 146 Gd and 118 Ba. 

Comparing current patterns shown in Fig. [51 one can 
conclude that for a system of non-interacting fermions, 
the total current, being the sum of the single-particle 
currents (see Eq. (fT7|l ). is, in general, quite complicated. 
This is a consequence of the fact that the localization, the 
strength and the structure of the current vortices created 
by a particle in a specific single-particle state depend on 
its nodal structure (see Ref. [43J and Sec. IIIC in Ref. 
Q). In this respect it is important to mention the results 
of Ref. [43| which showed that Coriolis induced current 
for a single-particle in a slowly rotating anisotropic har- 
monic oscillator potential has, in fact, a rather simple 
structure. It exhibits a number of localized circulations 
with precisely predictable centers and sense of rotation. 
The centers of the circulations are found at the nodes 
and peaks of the oscillator eigenfunctions, thus, forming 
a rectangular array somewhat similar to a crystal lattice. 

The wavefunction of the CRMF approach is more com- 
plicated than that of a rotating anisotropic harmonic os- 
cillator because of the presence of spin-orbit interaction 
and the split of the wavefunction into large and small 
components. Moreover, there are magnetic potential in- 
duced currents in addition to Coriolis induced ones. How- 
ever, the analysis of single-particle vortices in rotating 
nuclei in general confirms the observations made in Ref. 
[43| . The typical features of the single-particle currents 
in the CRMF approach are considered below on the ex- 
ample of three neutron single-particle states occupied in 
the yrast SD configuration of 152 Dy. The current and 
density distributions of these states are shown in Fig. 
O Let us first consider the i/[642]5/2+ state. The com- 
parison of Figs. IHt and[3Ji reveals that the rotation of a 
nucleus considerably increases the currents in this state. 
On the other hand, the density distribution is almost 
unaffected by rotation. The rotation of a nucleus also 
leads to a change in the structure of the circulations. At 
il x — 0.0 MeV, there are three weak circulations cen- 
tered around the nodes at (y = fm, z — fm) and 
(y = fm, z w ±4 fm); they are due to magnetic po- 
tential. Only two much stronger circulations are visible 
at Q x — 0.5 MeV: they are centered around the nodes 
located at (y = fm, z rs ±2.5 fm). This change of the 
structure of vortices can be attributed to additional cur- 
rents produced by the Coriolis term as well as to the 
change of the structure of wave function with increasing 



rotational frequency. The wave function in terms of two 
largest components has the 86%[642]5/2+5%[633]5/2 
and 63%[642]5/2+13%[651]3/2 structure 2 at Q x = 0.0 
MeV and fi x = 0.5 MeV, respectively. 

Even much large changes are induced by rota- 
tion into the structure of the ^[411]l/2 + state. 
The wave function in terms of two largest com- 
ponents has the 57%[411]l/2+23%[651]l/2 and 
84%[411]l/2+13%[411]3/2 structure at Q x = 0.0 
MeV and fl x — 0.5 MeV, respectively. One can see 
that the AN = 2 interaction, leading to a considerable 
admixture of the [651] 1/2 component into the structure 
of wave function, plays very important role at no 
rotation. The change in the wave function induced by 
rotation leads to a considerable changes both in the 
nodal structure of density distribution and in the current 
distribution (compare Fig. with Fig.[3f). 

The wave function of the z/[770]l/2 _ state is changed 
considerably by the rotation: its structure in terms of two 
largest components is 62%[770]l/2+17%[761]l/2 at ti x = 

0. MeV and 39%[770]l/2+28%[761]3/2 and Q x = 0.5 
MeV. The increase of rotational frequency does not lead 
to appreciable modifications in the density distribution 
but considerably decreases the strength of the currents 
and changes the shape of the circulations (see Figs.[S^,b). 
The latter is a consequence of additional Coriolis induced 
currents. It is interesting that for this state the currents 
show maximum strength at the densities far below the 
maximum densities. This most likely explains relative 
weakness of the currents in this state as compared with 
those in the ^[411]l/2 + state. On the contrary, for many 
single-particle states the strongest currents are seen at or 
close to local increases in the densities (see Figs. |Ht,d,f,g 
in the current manuscript and Fig. 8 in Ref. [5]). 

Our calculations show that the moments of inertia of 
the SD and HD configurations are very close to rigid- 
body values (Sec. IIV|) . However, the intrinsic currents 
show the dramatic deviations from rigid rotation. Usu- 
ally the deviations from the rigid-body moment of inertia 
imply that the flow pattern must substantially deviate 
from the current of a rigidly rotating mass distribution, 

1. e. there are strong net currents in the body-fixed frame 
[32| . However, the opposite is not true: the closeness of 
the moments of inertia to rigid body value does not nec- 
essary implies that the current distribution should corre- 
spond to rigid rotation. On a microscopic level, the build- 
ing blocks of the total current, namely, the single-particle 
currents certainly do not have a rigid-flow character; on 
the contrary, they have the vortex-flow character (see Fig. 

IS). 

Earlier non-relativistic studies also point to above dis- 
cussed relations between currents and rigid body mo- 



2 The percentages show the weights of respective components of 
the wave function in the total structure of the wave function. 
Note that only two largest components of the wave function are 
displayed. 
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merits of inertia. For example, it was shown in Ref. [42j 
for Schrodinger equation that single-valuedness require- 
ment for the wavefunction implies non-existence of rigid- 
flow in a quantum fluid. Furthermore, it was demon- 
strated for a system of independent particles employing 
cranked harmonic oscillator potential that the current is 
not of the rigid-flow type even when the moment of iner- 
tia assumes the rigid-body value ([I2> see also Ref. [4ljj ~) . 

Current distributions shown in Figs. [8] and [9] are typical 
for collective rotation around the a;— axis perpendicular 
to the symmetry axis. Note that the alignment of the an- 
gular momentum vector of a particle is specified along the 
x-axis in one-dimensional cranking approximation (see 
also discussion in Sec. IIIC of Ref. [5|). The = 1/2 
orbitals are aligned with the axis of rotation (x-axis) al- 
ready at no rotation. As a result, the single-particle an- 
gular momentum vector of the = 1/2 orbitals performs 
the precession around the x-axis, thus orienting the cur- 
rents predominantly in the y — z plane. In addition, the 
= 1/2 orbitals show vortices which are concentrated in 
the central region of nucleus. For the configurations with 
O 7^ 1/2, this mechanism of alignment becomes active 
only when the rotation sets up. Moreover, with increas- 
ing O, the densities and currents are pushed away from 
the axis of symmetry of the nucleus toward the surface 
area (Figs. GJfc,d and Fig. 8 in Ref. @). 



VI. FREQUENCY AND CONFIGURATION 
DEPENDENCES OF THE IMPACT OF 
NUCLEAR MAGNETISM ON THE MOMENTS 
OF INERTIA 

In this section, the frequency dependence of the impact 
of NM on the moments of inertia is studied using consid- 
erable number of SD and highly-deformed configurations 
in 60 Zn obtained in unpaired CRMF calculations. The 
properties of yrast SD band in this nucleus were well de- 
scribed in this formalism above band crossing which takes 
place at O^ ~ 1 MeV [HHfil, while the CRHB+LN for- 
malism gave good description of this band in the band 
crossing region |3l| . The neutron routhian diagram for 
this configuration obtained in the calculations with the 
NLSH parametrization is shown in Fig. 1 of Ref. [ll|; the 
results with the NL1 parametrization are similar to the 
ones obtained with NLSH. All proton and neutron states 
below the Z = 30 and N = 30 SD shell gaps are occupied 
in this configuration (note that proton routhian diagram 
is similar to the neutron one). The configurations are 
labelled by the shorthand notation [n,p], where n (p) is 
the number of occupied gg/ 2 neutrons (protons). In this 
notation, the yrast SD band has the [2,2] configuration. 
Excited configurations under consideration are built by 
means of proton or/and neutron particle- hole excitations 
across the Z = 30 and N = 30 SD shell gaps. 

The results of calculations for contributions of NM into 
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(2) 

dynamic (A J NM _ contr ) and kinematic 
moments of inertia are shown in Fig. 1101 At low fre- 



— contr 



10: (Color online) The contributions of NM to the dy- 
namic (J' 2 ') (panel (b)) and kinematic (J^) (panel (a)) mo- 
ments of inertia as a function of rotational frequency. Differ- 
ent color/line types are used for different groups of configu- 
rations characterized by the occupation of g g / 2 protons and 
neutrons. 



quencies, the average contribution of NM into kinematic 
moment of inertia is slightly larger than 20% (Fig. fTOk ) 
and the Ajj^jL contr quantities show considerable de- 
pendence on configuration. The origin of the latter ob- 
servation can be traced back to the specific features of 
some occupied single-particle orbitals. Let us consider as 
an example the [2,2] configuration. At low frequencies, 

the A J^m -contr values f° r this configuration are consid- 
erable higher than the A J^m -contr values averaged over 
all calculated configurations. This is due to the fact that 
upsloping branches of the proton and neutron [440]l/2 + 
orbitals (in the O^ = 0.0 — 0.7 MeV range, see Fig. 1 
in Ref. [11(), characterized by the expectation values of 
the single-particle angular momentum (j x )i strongly af- 
fected by NM, are occupied at O^ < 0.6 MeV. At frequen- 
cies O^ ~ 0.8 MeV, these orbitals strongly interact with 
proton and neutron [431]3/2 + orbitals and exchange the 
character of the wavefunction. This leads to unpaired 
band crossing (see Ref. |ll|) which is seen in consider- 
able changes of A J^_ contr and A J^ ] M _ contr quantities. 
The band crossing process is completed above O^ = 1.1 
MeV, where the orbital labeled as [440]l/2 + is downslop- 
ing as a function of rotational frequency (see Fig. 1 in Ref. 

[111]). At these frequencies, the AJ^| / _ contr quantity 
for the [2,2] configuration is slightly below the value of 

^NM-contr averaged over all calculated configurations 
(Fig. [TUk). Note that this unpaired band crossing is not 
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FIG. 11: (Color online) The contributions of NM to /32-deformation (panel (a)), 7-deformation (panel (b)), mass hexadecapole 
moment Q40 and total energy E (panel (d)) as a function of angular momentum 7. Note that the results on panels (b) and (d) 
are shown in absolute values. 



active in the [1,1] configurations because neither proton 
nor neutron [440]l/2 + orbitals are occupied. The cal- 
culations also suggest that it is considerably suppressed 
in the [3,3] configurations due to the changes in the de- 
formations and currents induced by the occupation of 
third (79/2 orbital both in proton and neutron subsys- 
tems. However, the presence of this crossing is still visi- 

( 2^ 

ble (especially, in the AJ y N ' M _ contr quantity) in some [2,3] 
configurations. 

With increasing rotational frequency, the average con- 
tribution of NM into kinematic moments of inertia de- 
creases and it falls below 15% at £l x ~ 2.5 MeV (Fig. 
ITUa). In addition, the configuration dependence of the 

^■■^NM-contr quantities is weaker than the one at low fre- 
quencies. At these frequencies, the majority of occupied 
single-particle orbitals are either completely aligned or 
very close to complete alignment. However, NM do not 
modify the expectation values of the single-particle angu- 
lar momenta (j x )i of completely aligned orbitals [lij]. As 
a result, only remaining orbitals, which are still aligning, 
contribute into A J^M-contr- The combined contribution 

of these orbitals into A J^M-contr * s smai l er than the one 
at lower frequencies because the alignment of these or- 
bitals is not far away from complete. 

The impact of NM on the dynamic moments of inertia 
is shown in Fig. 110b and it clearly displays much more 
complicated pattern as compared with the impact of NM 



on the kinematic moments of inertia. The irregularities 

(2) 

in the AJ]y M _ contr quantities are related to the band 

(2) 

crossings. For example, the dip in the A J^M-contr va ^" 
ues of the [2, 2] configuration at ft x ~ 0.9 MeV is caused 
by the unpaired band crossings which take place at differ- 
ent frequencies in the calculations with and without NM. 
Similar deviations from smooth trend as a function of 
rotational frequency are visible in other configurations. 
However, one can see that for some configurations the 
contribution of NM into dynamic moment of inertia is 
a smooth function of rotational frequency over extended 
frequency range. In this frequency range, the configura- 
tions remain unchanged. It is interesting that for some 
of these configurations the contributions of NM into dy- 
namic moments of inertia are either close to zero or even 
negative; such features have not been seen in the previ- 
ous analyzes of the impact on time-odd mean fields on 
the dynamic moments of inertia [H, 0, H, H| . 

The impact of NM on other physical observables of in- 
terest is shown in Fig. QTJ in which the results of the NM 
and WNM calculations are compared as a function of to- 
tal angular momentum. One can see that the quadrupole 
deformations P2 (Fig. ITTa ) obtained in the calculations 
with and without NM differ by less than 3%. The only 
exception is configuration B for which this difference 
reaches 7%. The difference in mass hexadecapole mo- 
ments Q40 obtained in the calculations with and without 
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NM is larger but typically below 10% fFig.fTTb'k the only 
exception is the configuration B for which this difference 
reaches 20% at I ~ 20H. The 7-deformations obtained in 
the calculations with and without NM differ by less than 
1.5° (Fig. [Tib). The only significant difference is seen in 
the total binding energies (Fig.lTTH). where the NM solu- 
tion is more bound than the WNM solution. This effect, 
which is due to the modifications in the moments of iner- 
tia induced by NM, is very large: additional binding due 
to NM reaches 7-8 MeV at spin / = 30h. These system- 
atic results are consistent with the ones obtained in the 
previous studies of single SD configuration in 152 Dy Q 
and single terminating configuration in 20 Ne [lij]. They 
also give a hint why the cranked models based on the 
phenomenological potentials like Woods-Saxon or Nils- 
son, which do not include time-odd mean fields [3J , are so 
successful in the description of experimental data. When 
considered as a function of spin the deformation proper- 
ties of the rotating system are only weakly affected by 
time-odd mean fields, and the proper renormalization of 
the moments of inertia [26| takes care of the E versus 
angular momentum curve. 



VII. PARAMETRIZATION DEPENDENCE OF 
THE CONTRIBUTIONS OF NM TO THE 
MOMENTS OF INERTIA 

It was shown in Ref. [5] that additional binding due 
to NM in one-particle states only weakly depends on the 
RMF parametrization; this is also seen in the analysis 
of terminating states in Ref. [lij]. In this context, it is 
important to understand how the contributions of NM to 
the kinematic and dynamic moment of inertia depend on 
the RMF parametrization. 

The dependence of the dynamic moments of inertia on 
the RMF parametrization has earlier been analyzed on 
the example of the SD bands in 151 Tb and 143 Eu in Ref. 
[l(| and in 58 Cu and 60 Zn in Ref. [n| employing the NL1, 
NL3 [13] and NLSH [H[ parametrizations of the RMF 
Lagrangian. The latter study includes also the results of 
calculations for kinematic moments of inertia. Additional 
calculations for these nuclei have also been performed 
with NL3* [33[ and NLZ [H] parametrizations for the 
current manuscript. As follows from these results, the 
kinematic and dynamic moments of inertia depend only 
weakly on the parametrization of the RMF Lagrangian. 
Indeed, at given frequency all the results for kinematic 
[dynamic] moments of inertia fit into the window which 
have a width equal to approximately 5% (ss 8%) [approx- 
imately 6% (« 10%)] of the value of kinematic [dynamic] 
moment of inertia in the A ~ 150 (^4 ~ 60) region of 
superdeformation. The larger spread of calculated val- 
ues in the A ~ 60 mass region are most likely due to (i) 
larger softness of potential energy surfaces in these nu- 
clei as compared with the ones in the A ~ 150 region of 
superdeformation and (ii) to larger relative importance 
of each particle and, thus, model uncertainties in the de- 



scription of their single-particle energies. 

Fig. [12] shows the dependence of the contributions of 
NM to the kinematic and dynamic moments of inertia 
on the parametrization of the RMF Lagrangian. For 
simplicity of comparison, these quantities are normal- 
ized to those obtained in the calculations with the NL1 
parametrization. Very weak dependence (within 5% win- 
dow with respect of the NL1 results) of the contribu- 
tion of NM to the kinematic moment of inertia on the 
RMF parametrization is seen in whole frequency range 
in 152 Dy (Fig.[JJb) and at frequencies Q x > 0.75 MeV in 
60 Zn (Fig. [T2"k). In the latter nucleus the deviation from 
the NL1 results reaches 10% at lower frequencies. The 
possible reasons for the larger dependence of calculated 
quantities on the parametrization in 60 Zn has been dis- 
cussed above. On the other hand, the deviations from 
the NL1 results are larger for the dynamic moments of 
inertia. These deviations can be as large as 8% at highest 
frequencies in the yrast SD configuration of 152 Dy (Fig. 
IT2H) and as large as 20% in the yrast SD configuration 
in 60 Zn (Fig. \TIb). Considering that the dynamic mo- 
ment of inertia is related to the second derivative of the 
total energy with respect of spin, a larger dependence of 
the dynamic moment of inertia on the parametrization is 
expected. 

These values can be used to estimate the uncertainty 
in the definition of the moments of inertia in the CRMF 
calculations due to the uncertainty in NM. The latter is 
related to the dependence of the [O nm - WNM ] param 
quantities (Fig. [T2")) on the RMF parametrization dis- 
cussed above. Dependent on nuclear system and con- 
figuration, the NM contribution to the total kinematic 
and dynamic moments of inertia is approximately 10-25% 
fSecs.lTVland lVip . Thus, the uncertainty of the definition 
of the absolute value of the total dynamic and kinematic 
moments of inertia due to the uncertainty in the defi- 
nition of NM is modest, being in range 0.5-5.0%. The 
fact that the moments of inertia of rotational bands of 
different structure in unpaired regime are well (typically 
within 5% of experimental data pi Iq|— Tl2l. l50l|) described 
in the CRMF calculations strongly suggests that NM and 
its impact on the moments of inertia is reasonably well 
defined in the CDFT theory. 



VIII. TERMINATING STATES 

The majority of rotational bands which do not have 
large deformation at spin zero will terminate in a non- 
collective terminating state at I max HE Hll IH3 3 - The 
regime of nuclear motion in terminatin g st ate is usually 
referred as 'non-collective rotation' [26Tl54|. This is be- 
cause of the fact that for an axially symmetric potential, 



Only recently the evidences for non-termination of some rota- 
tional bands at I m ax have been found |53|1 . 
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FIG. 12: (Color online) The contribution of NM to the physical observable O for specific parametrization of the RMF Lagrangian 
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quantity) normalized to the one obtained in the NL1 parametrization. The results of calculations 



are shown for yrast SD configurations in the Zn and Dy nuclei. Note that band crossing region is excluded in panel (b). 




FIG. 13: Neutron current distributions j n (r) in the intrinsic frame in the y — z plane (at x = 0.416 fm) (left panel), in the 
x — z plane (at y = 0.416 fm) (middle panel), and in the x — y plane (at z = 0.457 fm). These distributions are shown for the 
state of 47 V terminating at I — 17.5 + . The currents are plotted at arbitrary units for better visualization. The shape and size 
of the nucleus are indicated by density lines which are plotted in the range 0.01 - 0.06 fm -3 in step of 0.01 fm" 3 . Vertical and 
horizontal lines on left panel show the cross-sections at which the currents in the z — x plane (middle panel) and in the y — x 
plane (right panel) are plotted, respectively. 



the nucleon orbitals are not influenced by the rotation 
around the symmetry axis of this potential; thus, col- 
lective rotation about that axis is not possible. Non- 
collective rotation is also realized in the aligned states 
such as 'yrast traps' (or 'yrast isomers') |55l - [57j . The 
study of terminating states in the context of understand- 
ing of time-odd mean fields is of considerable interest 



because of several reasons. First, time-odd mean fields 
provide an additional binding to the energies of the spe- 
cific configuration, and this additional binding increases 
with spin and has its maximum exactly at the termi- 
nating state p3 |. This suggests that the terminating 
states can be an interesting probe of time-odd mean 
fields [lil l27l |58| provided that other effects can be re- 
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FIG. 14: Neutron current distributions j n (r) in the intrinsic frame in the y — z plane (at x — 0.407 fm). These distributions are 
shown for the state of 47 V terminating at I = 17.5 + . Left, middle and right panels show total, Coriolis induced and magnetic 
potential induced currents, respectively. See caption of Fig. [T3] for other details. 



liably isolated TlA l. Second, at the band termination, the 
NM does not modify either total angular momentum or 
the expectation values of the single-particle angular mo- 
menta (j x )i [3- Third, terminating state is a (multi)- 
particle+(multi)-hole non-collective state in which the 
angular momenta of all particles and holes outside the 
core are aligned along the symmetry axis. 

We will consider in this section the 7r(d3/2)i~5(/7/2)8 ® 
j/(/7/ 2 )| terminating state in 47 V, which has I m ax = 
17. 5 + , as an example. The structure of this state is 
given with respect of the 40 Ca core. This state is char- 
acterized by the largest impact of NM on the binding 
energies amongst terminating states studied in Ref. flij . 
It is nearly spherical with the quadrupole deformation 
/3 2 ~ 0.03 (Fig. 6 in Ref. QJ]). Our goal is to under- 
stand the impact of NM on the current distribution and 
microscopic origin of additional binding due to NM. 

In terminating states, the angular momenta of va- 
lence particles and holes are aligned along the symmetry 
axis (x-axis). As a consequence they perform precession 
around this axis, generating azimuthal currents with re- 
spect to the symmetry axis. This is illustrated in Fig. [TBI 
One can see two azimuthal circulations in the y — z plane: 
the circulation in the central region of nucleus is directed 
counterclockwise while the one in the surface region is di- 
rected clockwise. Fig. [TJ] shows total [left panel], Coriolis 
induced [middle panel] and magnetic potential induced 
[right panel] currents. One can see that surface circula- 
tion is generated by the Coriolis term, while the central 
circulation by the magnetic potential. The currents in 
the x — z and x — y planes are perpendicular to the x- 
axis (Fig. [TB")) . This clearly shows that the currents are 
azimuthal. 

In Ref. Q, the polarization effects induced by NM 
have been investigated in one- and two-particle configura- 
tions of odd and odd-odd non-rotating nuclei. Terminat- 
ing states differ significantly from these configurations. 
First, they are multi-particle+multi-hole configurations. 
For example, the 7r(d 3/2 )7^(/ 7/2 ) 4 (g)i/(/ 7 / 2 )f terminating 
state in 47 V has 8 particles and 1 hole outside the 40 Ca 



core. Second, the alignment of the angular momenta of 
these particles and holes generates considerable total an- 
gular momentum (/ = 17.5ft in the discussed terminating 
state of 47 V) aligned along the axis of symmetry; this mo- 
mentum is much larger than the one in odd and odd-odd 
nuclei studied in Ref. [f|. Third, terminating states are 
characterized by the azimuthal currents with respect to 
the symmetry axis, while the states in non-rotating nu- 
clei are characterized by the currents shown in Figs. 7-9 
of Ref. @. Thus, it is interesting to see how these differ- 
ences affect the polarization effects induced by NM and 
whether these polarization effects arc similar in nature 
for these two classes of non-collective states, namely, low- 
spin one- and two-particle configurations of non-rotating 
nuclei and high-spin terminating states. 

In order to facilitate the discussion, we split the total 
energy of the system (Refs. [6j,[9() into different terms as 4 



E,, 



part 



+ Ecor + E c 



-E 



TL 
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SL 
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Ec — E c 

- Ecoul, 



NL 



E 



TL 



(22) 



where E partl Ecor, and E cm represent the contributions 
from fermionic degrees of freedom, whereas the other 
terms are related to mesonic (bosonic) degrees of free- 
dom. In Eq. 



E 



part 



A 



(23) 



is the energy of the particles moving in the field created 
by the mesons (e^ is the energy of i-th particle and the 
sum runs over all occupied proton and neutron states), 



E, 



Cor 



i\jx\i) 



(24) 



4 We follow Refs. in the selection of the signs of the energy 

terms. 
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is the energy of the Coriolis term, 

E a = l -g a J d\a{r)[pP(r) +P :{r)] (25) 

is the linear contribution to the energy of isoscalar-scalar 
cr-field, 



EaNL = ^ / "' ' 



(26) 



is the non-linear contribution to the energy of isoscalar- 
scalar cr-field, 

El L = \a»J d 3 ru> (r) [pg(r) + p n v {r)] (27) 

is the energy of the time-like component of isoscalar- 
vector w-field, 



1 



E^ = -g p rf»rpo(r)[tf(r)-p?(r)] 



(28) 



is the energy of the time-like component of isovector- 
vector p-field, 



E 



SL 



1 

- 2 9u 



d 3 ru(r)[f(r)+j n (r)} (29) 



is the energy of the space-like component of isoscalar- 
vector w-field, 



E S p L = -\g p J d 3 rp(r)\j n (r)-f(r)} 



(30) 



is the energy of the space-like component of isovector- 
vector p-field, 



Ecoul = \ej d 3 rA (r)pP(r) 



is the Coulomb energy, and 



41A" 1/3 MeV 



(31) 



(32) 



is the correction for the spurious center-of-mass motion 
approximated by its value in a non-relativistic harmonic 
oscillator potential. 

Polarization effects induced by NM are investigated by 
considering NM impact on different terms of the total en- 
ergy (Eq. (|2"2l ). The results of this study are shown in 
Table H Similar to Ref. @, the Ej L and E p L terms are 
only weakly influenced by NM, and thus, they will not 
be discussed in detail. Somewhat stronger impact of NM 
is seen in the Ecoul, E a NL, and E^ L terms. Note that 
only last term was appreciably affected in low-spin con- 
figurations of odd and odd-odd nuclei in Ref. |5J . Much 
larger polarization effects are seen in the E part , E a and 
E^ L terms. The E a and E^ L terms depend only indi- 
rectly on time-odd mean fields through the polarizations 
of time-even mean fields induced by NM [5j . One should 



TABLE I: Impact of NM on different terms of the total en- 
ergy (Eq. IJ22J) m the state of 47 V terminating at I = 17. 5 + . 
Column (2) shows the absolute energies [in MeV] of different 
energy terms in the case when NM is neglected. Columns (3) 
and (4) list the changes AEi = E? M -Ef NM [in MeV] in the 
energies of these terms induced by NM in self-consistent [col- 
umn (3)] and perturbative [column (4)] calculations. Note 
that only nonzero quantities are listed in column (4). The 
results of calculations are obtained at rotational frequency 
£l x — 2.4 MeV. The energies of the Coriolis term Ecor and 
particle energy E par t depend on frequency. The latter takes 
place through the modifications of the energies of single- 
particle states with frequency (see, for example, Sec. 3.8 and 
Fig. 6a in Ref. |2j6j for more details on a construction of termi- 
nating state) . However, the sum Ecor + E part do not depend 
on frequency. Other remaining quantities shown in column 
(1) and the AEi quantities shown in columns (3) and (4) are 
frequency independent. 



Quantity 


j^WNM 


AEi 


AEf ert 


1 


2 


3 


4 




-1217.668 


-15.5 


-7.296 


Ecor 


42.0 







E a 


-6381.875 


-85.231 




EaNL 


109.368 


-2.815 






5371.721 


79.291 




e5 l 


0.0 


-3.51 


-3.51 


zpTL 

E P 


0.549 


-0.002 




E P 


0.0 


-0.038 


-0.038 


Ecoui 


102.168 


0.515 


0.240 


Ecm 


-8.521 


0.0 




Etot 


-386.121 


-3.704 


-3.983 



keep in mind that only the E a + E^ L quantity has a deep 
physical meaning, as it defines a nucleonic potential; this 
sum is modified by NM on -5.9 MeV. 

Comparing these results with those presented in Ref. 
Q , one can conclude that polarization effects for different 
total energy terms in terminating state under study are 
stronger by at least one order of magnitude than in low- 
spin one- and two-particle configurations of non-rotating 
nuclei. This is a consequence of the fact that all parti- 
cles (8) and holes (1) outside the 40 Ca core participate 
in building the total angular momentum and currents in 
terminating state, while only one (two) particle(s) par- 
ticipate in generating the currents in non-rotating odd 
(odd-odd) nuclei 5]. Despite that the relative impact of 
NM on different terms of the total energy is, in general, 
similar in these two classes of non-collective states (com- 
pare Table I in the present manuscript with Tables I, II, 
and IV in Ref. Q). 

Total modifications of the energies due to NM in the 
mesonic sector are -11.79 MeV. Only one third of these 
modifications comes from the terms (E^ L , E^ L ) which 
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directly depend on nucleonic currents, whereas the rest 
from the modifications of time-even mean fields induced 
by NM. 

It is interesting to compare the results of self-consistent 
and perturbative calculations 5 . The AEf ert = Ef M - 
fiWNM q uan tities will be used for simplicity in further 
discussion. These quantities are shown in columns 3 and 
4 of TableU The AE a , AE aNL , AE^ L and AEj L quan- 
tities are zero in perturbative calculations because time- 
even fields are fixed in these calculations. The AE^ L and 
AE s p L are the same in self-consistent (column 3) and per- 
turbative (column 4) calculations because the E^ L and 
Ej L terms depend only on time-odd mean fields, which 
are the same in the parts of the calculations that include 
NM. Particle energies E part are strongly modified by 
NM in self-consistent calculations; they change by —15.5 
MeV. Perturbative calculations show that only one half 
of AEp ar t is coming directly from time-odd mean fields 
(see Sec. IVB of Ref. [5j for more details on this mecha- 
nism) , the rest is due to polarization effects in time-even 
fields induced by NM. The same is true for the Coulomb 
energy term E Co ui- 

It is evident from Table Q] that 

AE self -const ^ AE P£t_ ^ 

Note that the superscript 'self-const' and 'pert' refers to 
fully self-consistent and perturbative results. The anal- 
ysis of polarization effects in other terminating states 
of the A ~ 40 mass region shows the same relation. 
These results clearly indicate that the additional bind- 
ing due NM (the E NM - E WNM quantity) is defined 
mainly by time- odd fields and that the polarization effects 
in fermionic and mesonic sectors of the model cancel each 
other to a large degree. The same result has been earlier 
obtained in the analysis of non-rotating nuclei in Ref. Q . 



IX. SIGNATURE-SEPARATED 
CONFIGURATIONS 

Signature separation phenomenon induced by time- 
odd mean fields has been found earlier in excited 4- 
particle SD configurations of 32 S [HI, Hfll, and very re- 
cently in 2-particle configurations of nonrotating odd-odd 



5 Fully self-consistent calculations with NM provide a starting 
point for perturbative calculations. Using their fields as input 
fields, only one iteration is performed in the calculations with- 
out NM: this provides perturbative results. Time-even mean 
fields are the same in both (fully self-consistent and perturba- 
tive) calculations. Then, the impact of time-odd mean fields on 
calculated quantities (for example, different terms in the total 
energy (Eq. 122> 1 is defined as the difference between the values 
of this quantity obtained in these two calculations. In this way, 
the pure effects of time-odd mean fields in fermionic and mesonic 
channels of the model are isolated because no polarization effects 
are introduced into time-even mean fields [a|. 



nuclei in Ref. Q ■ It reveals itself in a considerable energy 
splitting of the r tot = +1 and r tot = — 1 branches of the 
configurations which have the same structure in terms 
of occupation of single-particle states with given Nilsson 
labels. Such a signature separation could not have been 
obtained in phenomenological cranking models, such as 
the ones using the Woods-Saxon or Nilsson potentials, 
since time-odd mean fields are absent in these models. 

However, the description of rotating N w Z nuclei re- 
quires isospin projection [6l| which can modify above 
mentioned results. Since this projection is beyond the 
current framework, we concentrate at the nuclei away 
from the N = Z line. The analysis of Ref. [f| shows 
that signature separation is expected also in such nuclei, 
but it is weaker as compared with the one seen in the 
nuclei around the N = Z line. Unfortunately, the sur- 
vey of odd-odd A — 20 — 52 nuclei (some of which were 
studied in Ref. |f|) does not reveal experimental bands in 
the nuclei away from the N = Z line in which signature 
separation is expected. 

Fig. [T5] shows that signature separation phenomenon 
can also be present in heavier nuclei. This figure shows 
the results of calculations for odd-odd Eu isotopes in 
which odd proton occupies fixed 7r[532]5/2 + state, and 
odd neutron occupies different neutron states of the 
r = ±i signatures along the isotope chain. Additional 
binding due to NM (the E NM -E WNM quantity, see Ref. 
5] for more details) is shown for total proton-neutron 
configurations with different total signatures r to t = in- 
significant signature separation (on the level of 100- 
150 keV) is seen in the tt[532]5/2 eg) z/[523]5/2 ( 156 Eu), 
?r[532]5/2® z/[642]5/2 ( 158 < 160 Eu), tt[532]5/2 ® i/[633]5/2 
( 196 Eu), and tt[523]5/2 ® z/[752]5/2 ( 204 Eu) configura- 
tions. Either r = — lorr = +l states can be more 
bound in signature separated configurations of Eu iso- 
topes (Fig. fT5|) . This depends on mutual orientation of 
proton and neutron currents induced by odd proton and 
odd neutron; the state with the same orientation of these 
currents is more bound. 

Fig. [16] illustrates that four rotational sequences (two 
with total signature r to t = +1 and two with r to t = —1) 
can be built in the 2-particle configurations ir\a > (r = 
±z) ® v\b > (r = ±i) (where \a > and \b > indicate the 
blocked proton and neutron Nilsson states, respectively) 
of odd-odd nuclei. For the case of 158 Eu we consider 
2-particle configurations based on the \a > = [532]5/2 
and \b >= [642]5/2 states. In the WNM calculations, 
the r to t — +1 and r to t = — 1 configurations are almost 
degenerate in energy up to spin / ~ 10h (Fig. [TB|) . On 
the contrary, there is a considerable signature separation 
(« 150 keV) due to time-odd mean fields between these 
configurations in the calculations with NM. This feature 
is a strong spectroscopic fingerprint of the presence of 
time-odd mean fields. Note that rotational sequences A 
and B undergo unpaired band crossings at / ~ 20H. 

Unfortunately, the experimental data on odd-odd Eu 
nuclei also do not reveal the configurations discussed 
above in which the signature separation is expected. This 
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FIG. 15: The impact of NM on binding energies of the con- 
figurations in odd-odd Eu (Z = 63) nuclei at S} x = 0.0 MeV. 
The E NM - E WNM quantity is shown for different total sig- 
natures r. Proton configuration has blocked 7r[532]5/2 + state 
for all isotopes, while the neutron configuration changes as a 
function of mass number. Configurations are labeled with the 
Nilsson labels of blocked odd neutron state; configurations at 
and to the right of the Nilsson label up to the next Nilsson 
label have the same blocked neutron state. 

situation may be resolved by a systematic search of sig- 
nature separated configurations both in the experimen- 
tal data on odd-odd rare-earth nuclei and in the model 
calculations. The best way to confirm the existence of 
this phenomenon would be to find in model calculations 
the configurations of odd-odd nuclei which show no sig- 
nature splitting in the absence of time-odd mean fields 
and measurable signature separation in the presence of 
time-odd mean fields, and then to find their experimen- 
tal counterparts which show signature separation. How- 
ever, such an investigation is definitely beyond the scope 
of the current study. The difficulty of such a study is 
also underlined by the fact that existing interpretations 
of two-quasiparticle configurations in odd-odd nuclei are 
based on Woods-Saxon or Nilsson potentials. In these 
potentials, the signature degeneracy is considered to be 
a strong fingerprint of specific configurations. However, 
time-odd mean fields in EDF provide additional mecha- 
nism of breaking the signature degeneracy, so the exper- 
imental data on such configurations has to be reanalyzed 
in density functional calculations. 

Although the model calculations clearly indicate the 
important role of time-odd mean fields in creating signa- 
ture separation phenomenon in odd-odd nuclei, the direct 
comparison with experiment will be complicated by the 
number of model limitations which are related to 

• the presence of residual proton- neutron interaction 
of unpaired proton and neutron, 

• the coupling scheme of angular momenta vectors of 
unpaired proton and neutron at low spin. 

In odd-odd nuclei the angular momenta of unpaired 
proton and unpaired neutron in 2-quasiparticle config- 



urations can be coupled either in parallel or antiparallel 
fashion, namely into K > = Q p + £l n , and K < = |f2 p — fi n |, 
where fl p ( n ) represents the projection of single quasipar- 
ticle angular momentum of proton (neutron) on the axis 
of symmetry. For example, in 158 Eu this will lead to 
rotational sequences with K < = and K > = 5. The 
degeneracy of the bandheads of the Ky< doublet pair 
(called Gallagher- Moszkowski doublet [62|]) is lifted by 
inclusion of the residual proton-neutron interaction and 
also by the zero-point rotational energy. Relative en- 
ergy ordering of the K > and K < bands is determined 
by the empirical Gallagher-Moszkowski (GM) rule which 
places the spin-parallel band lower in energy than its 
spin-antiparallel counterpart 1631 in odd-odd nuclei (and 
vise versa in even-even nuclei [631 ] ) and has only few ex- 
ceptions [64l l65j . Another important consequence of the 
residual interaction of unpaired nucleons is the observed 
shift of the odd- and even-spin rotational levels relative 
to each other in the K = bands; this feature is generally 
referred to as the Newby or odd-even shift [6|| . 

Residual proton-neutron interaction of unpaired nu- 
cleons is neglected in the cranking models; we are not 
aware about any publication which includes it. So, nei- 
ther Gallagher-Moszkowski splittings nor Newby shifts 
can be described in the current calculations. It is also 
necessary to recognize that 2-quasiparticle configurations 
in odd-odd and even-even nuclei show a daunting com- 
plexity due to the high density of states and the large 
number of couplings and interactions possible. The prob- 
lem of the description of the Gallagher-Moszkowski split- 
tings and Newby shifts is far from being settled even 
in the framework of conventional particle+rotor model 
[6314651 l67l l68j. For example, the residual interaction of 
unpaired proton and neutron in odd-odd nuclei shows 
pronounced dependence on the mass region under study 
[6J, [65|. It is even more difficult to understand why 
in 2-quasiparticle configurations of the rare-earth region 
different residual interactions are required to describe 
the interaction between unpaired proton and neutron 
in odd-odd nuclei and between unpaired protons (neu- 
trons) in even-even nuclei [63| despite the expectations 
that they should be the same due to isospin symme- 
try. To our knowledge, self-consistent description of 
Gallagher-Moszkowski splittings and Newby shifts has 
been attempted only in the framework of the rotor+two- 
quasiparticle model based on Skyrmc Hartree+Fock ap- 
proach in Ref. [69j. 

At zero rotational frequency the angular momenta of 
odd proton and odd neutron are aligned (parallel or anti- 
parallel) with the symmetry axis which leads to band- 
head states with K > — fl p + fl n and K < = \£l p — f2 n |. 
However, in one-dimensional cranking approximation nu- 
clear configuration on top of which rotational sequence 
is built does not depend on coupling of Q p and £l n . This 
is well known (although seldom stressed) deficiency of 
one-dimensional cranking approximation. However, with 
increasing rotational frequency the angular momenta of 
odd proton and odd neutron start to align with the axis of 
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FIG. 16: (Color online) The energies of calculated bands in 
158 Eu based on the tt[532]5/2 ± ®^[642]5/2 ± two-particle con- 
figurations as a function of angular momentum with rigid ro- 
tor reference subtracted. The results of the calculations with 
and without NM are shown. 



rotation which is perpendicular to the axis of symmetry. 
Although it is tempting to employ tilted axis cranking 
(TAC) approximation for the description of the combina- 
tion of these two angular momenta coupling schemes at 
low spin, this does not resolve the problem of the descrip- 
tion of signature separation since signature is no longer 
good quantum number in the TAC approximation [70|. 
On the contrary, one-dimensional cranking approxima- 
tion used in the current manuscript has a clear advantage 
that it properly accounts for the alignments of valence 
particles and holes along the axis of rotation at medium 
and high spins where I > K [57[, and thus provides cor- 
rect description of signature separation at these spins. 



X. CONCLUSIONS 

Time-odd mean fields (nuclear magnetism) have been 
studied in rotating nuclei in a systematic way within the 
framework of CDFT. The main results can be summa- 
rized as follows. 

• NM can considerably modify the band crossing fea- 
tures (crossing frequencies, the properties of kine- 
matic and dynamic moments of inertia in the band 
crossing region). In the calculations without pair- 
ing, these modifications depend on the underly- 
ing changes in the single-particle properties such 
as alignments and energies induced by NM. These 
effects are also active in the calculations with pair- 
ing. In addition, in the calculations with pairing 
the gradual breaking of high- j pairs proceeds faster 
in the presence of NM, which is reflected in a faster- 
decrease of pairing with increasing fi x . Thus we 
can specify this effect as an anti-pairing effect in- 
duced by NM. 



• Outside the band crossing regions, the contribu- 
tion of NM to the kinematic and dynamic mo- 
ments of inertia only weakly depends on the RMF 
parametrization. 

• It is shown for the first time within the self- 
consistent approach that the moments of inertia 
of super- and hyperdeformed configurations in un- 
paired regime come very close to the rigid-body 
values. Despite that the presence of strong vor- 
tices demonstrates the dramatic deviation of the 
currents from rigid rotation. On the contrary, the 
moments of inertia of normal-deformed nuclei de- 
viate considerably from the rigid-body value in the 
calculations without pairing. 

• Complicated structure of the currents in the rotat- 
ing systems of independent fermions is the conse- 
quence of the fact that total current is the sum of 
the single-particle currents. The single-particle cur- 
rents show vortices (circulations), the strength and 
localization of which depend on the single-particle 
state. 

• Within specific configuration the impact of NM 
on the binding energies reaches its maximum at 
the terminating state [14j . Underlying microscopic 
mechanism for additional binding due to NM at 
such states has the same features as those seen 
in low-spin one- and two-particle configurations of 
odd and odd-odd nuclei [5(. However, the mag- 
nitude of the effects is significantly larger. The 
perturbative results clearly indicate that additional 
binding due NM at terminating states is defined 
mainly by time-odd fields and that the polariza- 
tion effects in fermionic and mesonic sectors of the 
model cancel each other to a large degree. 

• Signature-separation phenomenon in odd-odd nu- 
clei has been analysed in detail. It is shown that 
the effects neglected in the current approach such 
as the residual interaction of unpaired proton and 
neutron and the coupling scheme of angular mo- 
menta vectors of these particles at low spin con- 
siderably complicate quantitative description of the 
spectra of odd-odd nuclei. The best way to confirm 
the existence of this phenomenon would be to find 
(both in experiment and in calculations) the config- 
urations of odd-odd nuclei which show no signature 
splitting in the absence of time-odd mean fields and 
measurable signature separation in the presence of 
time-odd mean fields. 

Although time-odd mean fields affect different physical 
observables (see introduction in Ref. [5j for details), this 
investigation clearly shows that rotating nuclei still offer 
one of the best probes of this channel of density func- 
tional theories. This is because the impact of time-odd 
mean fields is significant representing on average 20% of 
kinematic and dynamic moments of inertia. In addition, 
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it shows appreciable variations with configuration, par- 
ticle number and rotational frequency; these variations 
provide a useful tool for a better test or definition of time- 
odd mean fields. Significant amount of the data on differ- 
ent types (normal- jl2| , superdeformed [l|, l^-flll l3ol [7l| , 
and smooth terminating [l|, [26|) of rotational bands in 
unpaired regime available in different mass regions offers 
a testing ground for time-odd mean fields. This data is 
also extremely useful for fitting the parameters of time- 
odd mean fields as needed, for example, in Skyrme energy 
density functionals, in which these fields are not well de- 
fined (Refs. d,[2!|). Our investigation, however, suggests 
that such fit has to be performed to a significant set of 



rotational structures representing different mass regions 
and different configurations and spanned over significant 
frequency range in order to minimize the dependence of 
the fit parameters on the choice of experimental data. 
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